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ABSTRACT 


Traction tests were performed on RP1, a common kerosene based rocket 
propellant. Traction data on this fluid are required for purposes of 
turbopump bearing design, using codes such as SHABERTH. To obtain the 
traction data, an existing twin disc machine was used, operating under the 
side slip mode and using elliptical contacts. The range of test variables 
were: contact peak Hertz stress from 1.0 to 2.0 GPa, disc surface speed from 

10 to 50 m/s, fluid inlet temperature from 30 to 70 °C, and with a contact 
aspect ratio of 1.7. 

The resulting traction curves were reduced to fundamental fluid property 
parameters using the Johnson and Tevaarwerk traction model. Theoretical 
traction predictions were performed by back substitution of the fundamental 
properties into the traction model. Comparison of the predicted with the 
measured curves gives a high degree of confidence in the correctness of the 
traction model. For purposes of input to the NASA SHABERTH program, the 
traction model was next used to predict the expected traction of RP1 under 
line contact conditions. 
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NOMENCLATURE 


Below follows a list of symbols used in the text and their units. 
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Fluid viscosity temperature parameter 

Fluid viscosity pressure parameter 
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Shear stress temperature constant 
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Pressure solidification temperature 

Viscosity solidification temperature 

Shear stress solidification temperature 
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Viscosity constant for non-linear thermal model 
Composite elastic modulus for the disc material 
Experimental traction viscosity constant 
Experimental time delay constant 
Dissipative function for traction model 
Dimensionless thermal resistance of the film 
Contact force in the x direction 
Contact force in the y direction 
Normal force on the contact 
Fluid shear modulus (uncorrected) 

Compliance corrected fluid shear modulus 
Johnson's elasticity parameter 
Shear modulus of the disc material 
Central film thickness in contact 
Dimensionless longitudinal slip variable 
Dimensionless longitudinal traction variable 
Elastic stress portion of J4 
Plastic stress portion of J4 
Thermal dimensionless longitudinal traction 
Contact aspect ratio (b/a) 

Thermal conductivity of fluid 
Thermal conductivity of disc material 
Fluid thermal conductivity constant 
Fluid thermal conductivity constant 
Calibration constant in side slip measurement 
Initial slope of fluid traction curve 
Initial slope of dry traction curve 
pressure 

Mean Hertz contact pressure 
Hertzian peak contact pressure 
Psuedo Peclet number 
Reduced Hertz pressure 
Kalker coefficient 

Equivalent radius of curvature for discs 
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i v 


s 


s 

t 

u 

AU 

V 0 

AV 

AX 

YI 


Auxiliary variable used in elastic/plastic model 

Roller slip 

Time 

Rolling speed 

Slip velocity in x direction 

Viscosity constant 

Side slip velocity of the discs 

Small displacement of the displacement transducer 

Inlet shear heating factor 


s 
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.Pa.s] 
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Cm] 

r l 


GREEK SYMBOLS 


P 

Q 

0 o 

$ 



■s 

T c 

T .o 

7 

/* 

V 

Vo 


Side slip angle 

Temperature of the fluid 

Shearplane temperature 

Inlet or blank temperature 

Time delay parameter 

Density of the disc material 

Hertzian contact shape parameters 

Fluid shear stress 

Non linear shear stress parameter 

Limiting shear strength of fluid 

Limiting shear strength at inlet temperature 

Shear strain rate 

Peak traction coefficient 

Viscosity of fluid in the contact 

Viscosity constant 
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°C 

°C 

°c 
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1-0 INTRODUCTION 


Traction or friction plays a major role in today's technological society 
in that it holds one of the keys to reduce our overall energy consumption, and 
thereby the dependence on unreliable sources of this energy. Friction and 
traction indicate the resistance to relative motion of two 'contacting' 
bodies. The term traction and friction have the same meaning in a 
tribological sense, however friction is used when this resistance is 
undesirable and traction is used when it is desirable. The mechanical 
components in which friction and traction are important are rolling element 
bearings, gears, cam and tappets, and traction drives. Because these devices 
almost always operate in a wet or fluid lubricated environment, the traction 
or friction is mostly governed by the particular fluid that is used. In the 
first three devices friction is the key source of inefficiency and because of 
the multitude of bearings and gears in service, a small reduction in these 
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losses can amount to phenomenal savings in energy. Rolling element bearings 
actually have rather interesting requirements for friction or traction in that 
at low to medium speeds friction should be low, but at high speeds traction 
should be high to ensure that the rolling elements operate at the correct 
velocities. 

In most mechanical devices where friction or traction plays an important 
role, the lubricating fluid is chosen to optimize the energy efficiency of the 
device, however, this is not possible in all cases. For example, the bearings 
used in rocket motor turbopumps are lubricated with the process liquid, may 
this be a mineral-based low-viscosity fuel such as RP1, or a cryogenic 
substance such as liquid oxygen or hydrogen. For applications such as these 
it is critically important that the rheological properties of the working 
fluid be known and modeled such that a good deal of confidence exists about 
using the data beyond the range of observation. The data thus obtained can 
then be used in bearing design codes such as SHABERTH to design for adequate 
bearing performance for a given mission. 

1-1 Prior Traction Investigations 

There has been a lot of activity in the area of traction research, both 
in the past and recently. Notable contributions have come from Clark et al 
[1951], Smith [1965], Smith et al [1973], Johnson and Cameron [1967], Niemann 
and Stoessel [1971], and more recently Johnson and Roberts [1974], and Johnson 
and Tevaarwerk [1977]. Some of these investigations were strictly 
experimental in nature and aimed at obtaining traction design data, while 
others were aimed at understanding the traction phenomena so that rheological 
models could be formulated. This latter research is of course ultimately 
aimed at relating fluid molecular properties to traction properties. Research 
by Johnson and Tevaarwerk [1977], Daniels [1978], Hirst and Moore [1980], and 
Alsaad et al [1978] is directed more towards this purpose. The reader is 
referred to an excellent review by Johnson, [1978] for further aspects of this 
topic. 

Current understanding of traction has led to traction models that 
describe the fluid shear behavior in terms of an elastic and a dissipative 
element. For purposes of mathematical tractability this dissipative element 
is taken to be plastic like in nature. This gives an adequate description of 
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the fluid behavior at conditions such as those encountered in rolling element 
bearings. This work has now been further expanded by developing a simple 
method to correct for thermal effects due to spin, Tevaarwerk [1981], and an 
overall thermal traction study, Tevaarwerk [1980]. 

As with all models, however, their usefulness is severely restricted if 
inadequate input data is available to the designer. This is especially so if 
fluids are used that are not normally desirable from a lubrication point of 
view, but are selected on the basis of other features. RP1 is such a fluid. 
Being a rocket propellant, based on a high viscosity kerosene, it does not 
have the traction data base that other more lubricous fluids have. Therefore, 
it was decided to test its traction behavior experimentally and to use 
existing rheological models to predict its traction under conditions different 
from those in the experiments. 

The investigation reported herein was performed by Battelle Memorial 
Institute, of Columbus, Ohio, under contract to the NASA Lewis Research Center 
through the Metals and Ceramics Information Center, a DOD "Information 
Analysis Center", Item 0001CJ, Control No. DLA 900-83-C-1744. The NASA 
technical project manager was C. Woods of the Launch Vehicle Technology Branch 
in the Space Propulsion Division. 
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2-0 EXPERIMENTS 

The various traction experiments were carried out on an existing twin 
disc test facility, shown in Figures 2-1, 2-2 and 2-3. The test facility was 
also used for the traction test results as reported in Tevaarwerk [1983] and 
[1985], Traction curves are obtained by using the side slip technique so 
that large traction transfer can be measured without the need for a large 
motor/generator set. It also has the advantage that the effects of bearing 
friction and rolling traction can be corrected for by the simple measurement 
of the longitudinal slip. This is particularly important when the expected 
levels of fluid friction are low, as is the case with the lubricant tested 
here. With longitudinal traction measurements the correction for bearing 
friction and rolling traction is often of the same magnitude as the measured 
traction with these lubricants. 

2-1 Description of Twin Disc Machine. 

For an extensive description of twin disc traction testers the reader is 
referred to the literature; Smith [1965] and Johnson and Roberts [1974]. 
Basically the machine consists of two discs, called the toroid (upper) and 
disc (lower). The lower disc is mounted in rolling element bearings through 
shafts and the only degree of freedom is one of rotation about this axis. The 
lower disc always has a transverse radius of curvature of infinity. The upper 
disc is contained in bearings that are mounted in the upper assembly. This 
upper assembly is suspended with elastic hinges such that only direct normal 
motion and axial motion is possible. The assembly will however always stay 
horizontal. The toroid (or upper) has curvatures to produce the desired 
contact geometry is arrived at. The upper assembly is constructed such that 
the toroidal axis can be tilted relative to the horizontal plane so as to 
introduce spin on the contact. It can also be skewed about the normal to the 
contact to introduce a side slip velocity. 

In order to achieve the range of pressures required for the traction 
data, a set of discs with a nominal aspect ratio of 1.7 is used. These discs 
are made of AISI-01 steel, hardened to 7.00 GPa, ground and polished to a 
surface finish of less than .05 pm RMS and with an out-of-roundness error of 
less than 5 pm. Between tests the discs are inspected for surface damage and, 
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if needed, reground and polished to bring them back up to specifications. The 
required normal load, obtained by a dead loading technique, can be calculated 
from the Hertz theory for elastic bodies in contact. The maximum contact 
normal stress P 0 is given by: 


( 2 - 1 ) 


where, 


F 

P 

E 


S 


R 


e 


*. x 


1 3 3F Z E' 2 

Wx K 2 


contact normal load, 

Hertzian contact stress, 
composite elastic modulus, 

(231 GPa for steel) 
equivalent disc radius, 

(1/R X + 1/R y ) _1 , and 

Hertzian contact shape factors. 


[N] 

[Pal 

[Pa] 

M 

[-1 


2-2 Instrumentation of the Traction Tester 

By suitably instrumenting the disc machine the relevant experimental 
parameters can be measured. In the experiments reported here the sideslip, 
sideways traction force, disc surface temperature, rolling velocity, and the 
amount of longitudinal slip are measured. The technique of measuring each of 
these variables will be discussed next. 

2-2-1 The Measurement of the Traction Force 

A ring dynamometer type load cell is used to measure the side slip force 
of the upper toroid assembly. The electrical signal from the load cell is 
conditioned for noise and amplified using common mode rejection techniques. 
The gain on the amplifier is adjusted so that a good range on the signal is 
measured for each test. Calibration of the load cell is done in situ by dead 
loading. This calibration is checked periodically. 
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2-2-2 The Measurement of Side Slip 

The skew angle is measured by using a direct current displacement 
transducer on the upper assembly and thereby measuring the rotation angle of 
this assembly. This skew angle gives the amount of side slip/roll ratio 
through the relationship: 

(2-2) Av/U = tan (/» , 


where, 

rad 
m/s 
m/s 


p = side slip angle, 

Av = side slip velocity, and 
U = rolling velocity. 


By measuring the amount of skew with the displacement transducer the side 
slip/roll ratio is obtained directly through: 

(2-3) Av/U = Kax . 

where Ax is the displacement of the transducer and K a scale factor. The 
electrical output of the displacement transducer is filtered in an R-C network 
to provide a low-pass signal. The maximum frequency response of the R-C 
network is 0.1 s. The displacement transducer is calibrated by rotating the 
upper assembly through a known angle and then calculating the amount of side 
slip for this angle. 

2-2-3 The Measurement of Rolling Velocity and Longitudinal Slip 

The rotational velocities of both the disc and the toroid are measured by 
using a MC 6840 frequency counter. The proximity probes for the control of 
the counters are mounted near the support bearings, see Fig. 2-3. Each shaft 
has a simple protrusion on it that produces a pulse once per revolution. This 
pulse is used to turn countdown counters on when controlled to do so. The 
next successive pulse from the proximity probe stops the countdown. When each 
of the two counters has completed its cycle, a flag is set and the register 
contents are read. From the difference between the original and the final 
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contents of the registers, the amount of time for a single revolution of each 
shaft is calculated. From this the angular speed for each shaft is 
determined, and hence the peripheral velocity of each roller can be calculated 
if the roller radius is known. The countdown frequency on the counters is 
selected such that an accuracy of at least 1: 10000 is obtained. By employing 
this method the slip of the rollers can be measured every five or six 
revolutions of the shafts. From the countdown and the undeformed roller radii 
the slip is calculated as follows: 

{ n ^ _ _ o U2 - U1 _ J R2 C2 W ( R2 A C2 ) r , 

(2_4) s ' 2 urn jt ' 2 V in ■ uJ / Vm rr J 1 J • 

This can be approximated if the amount of slip is small by the following: 

( 2 -s) s ■ v in' ' n ) • 

where , 

s = roller slip, 

U = rolling velocity, 

C = number of count down pulses, and 
R = roller radius. 

Suffix 1 denotes the toroid and suffix 2 denotes the disc. It should be 
stressed that the actual amount of slip will be somewhat different when two 
different roller materials are in contact because of the deformation of the 
rollers themselves. The errors thus introduced can be removed by measuring 
the rolling slip at very low rolling speeds. 

2-2-4 The Measurement of the Disc Surface Temperature 

In the analysis and reduction of the test data it is important that the 
disc temperature be used as the reference inlet temperature for the film 
thickness. This temperature can be measured by embedding a thermocouple 
directly below the surface of the disc and then to take this signal out 
through slip rings. This method is however not very practical when a number 
of different discs are involved and is also very costly from an installation 
point of view. With care the surface temperature can also be measured by 


[-] 

m/s] 

>3 
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using a trailing thermocouple that rides on the disc surface. The 
disadvantage of this technique is that it can be speed sensitive in its 
response because of frictional heating. 

The latter technique is employed here and care is taken to ensure that 
the contact force on the thermocouple is not excessive. A temperature 
reference bath junction ensures that the same reference level for the 
thermocouple is used at all times. The signal from the thermocouple is 
amplified using common mode rejection techniques to minimize the influence of 
electrical noise and other disturbances. Calibration is done by the boiling 
water method adjusted for sea level differences. This calibration is checked 
periodically. Only slight deviations are encountered. Because of the 
frictional heating at the junction/toroid interface a variation of about 2 C° 
is found in the signal between stationary discs and those rotating at a 
surface velocity of 120 m/s. The overall reproducibility of the temperature 
measurement is better than 2 C°. 

The temperature on the machine is regulated through the use of heaters 
and coolers on the test fluid. This test fluid is allowed to circulate freely 
before the start of a test series in order to bring the machine up to a 
uniform temperature. No specific effort is made however to maintain a given 
set point temperature during the test and the reason for this will become 
clear in the analysis of the results. 

2-3 Traction Measurements 

Traction curves are obtained by the slow rotation of the upper assembly 
from a positive value of side slip/roll ratio to a negative value. The 
signals from the above discussed transducers are fed into a digitizer from 
where they are led into a data logging computer for plotting and storage on 
magnetic media for future use. The computer automatically traces the force 
versus slip curve on the screen. By reversing the direction of rotation of 
the machine, a duplicate set of curves can be obtained. For each experiment 
500 data points are taken at fixed time periods of 0.2 seconds. Multiple data 
points are stored as separate entries. 

After the completion of a test series the data is recalled into memory of 
the computer and further manipulated. This manipulation consists of the 
averaging of the multiple entries, the filling in of any gaps in the data 
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through forward and backward interpolation, the comparison of the traces for 
the forward and reverse rolling direction, and the centering of the traces 
about the center lines. After the centering operation the data is smoothed by 
a ’ N ' point averaging technique for traction points after the peak traction 
points. For storage a geometric series is used so that the total traction 
trace is now represented by 40 data points for each measured variable. These 
traces are stored on magnetic media and used for further manipulation and data 
extraction at a later point. 

2-4 Traction Curves For RP1 


A survey of the literature revealed that very little traction data is 
available on low viscosity fluids such as RP1. This is especially so in the 
higher speed ranges of 10 m/s and above. Hence it was decided to perform 
traction tests on RP1 over a fairly broad speed, pressure, and temperature 
range. The nominal test matrix consisted of the following: 


Contact Pressure P 0 : 
Rolling speed U: 
Temperature 0 O : 

Slip AV/U: 
wl ab 

Spin — g — : 


1, 

1.3, 1.6, and 2.0, I 

;cpa; 

± 

10, ± 30, ± 

50, 

’m/s' 

25 

and 75, 


°C] 

-0 

.05 to 0.20, 

and 

[-] 

0. 


[-] 


This test matrix leads to 24 different traction test curves. The 
measured value being the traction coefficient F y /F z as a function of slip. 
The traction data measured on RP1 using the equipment and techniques as 
described above are shown in Appendix I. The lettered labels in each curve 
are used only for identification purposes and do not indicate an actual 
measured data point. 
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FIGURE 2-1. OVERVIEW OF HIGH PERFORMANCE SIDE SLIP TRACTION TESTER 
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FIGURE 2-2. RIGHT HAND VIEW OF HIGH PERFORMANCE SIDE SLIP TRACTION TESTER 
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3-0 THEORETICAL ANALYSIS 

In order to understand the required analysis of the experimental data it 
will be helpful to consider the following discussion of traction. The ability 
of a fluid film, trapped under high pressure in the elastically deformed 
region of two loaded curved elements, to transmit a tangential force from one 
element to the other is commonly referred to as friction or traction. The 
magnitude of this force depends on several variables such as :1) the contact 
kinematic conditions of slip, spin and sideslip, 2) the fluid present, 3) 
temperature, pressure and operating speeds. We will examine the traction 
behavior under simple slip only. 

Under conditions of increasing slip between the two elements, an 
increasing traction force is transmitted up to a certain limit at which point 
it will decrease with further slip. See Figure 3-1. 



FIGURE 3-1. TYPICAL TRACTION /SLIP CURVE 

There are three regions identified on this traction curve and the 
behavior in each of these regions can best be described by the Deborah 
number (De). For a simple Maxwell viscoelastic model this number is the ratio 
of the relaxation time and the mean transit time, see Johnson and Tevaarwerk, 
[1977]. 

(A) The linear low slip region. Thought to be isothermal in nature, it is 
caused by the shearing of a linear viscous fluid (low De) or that of a linear 
elastic solid (high De) . 

(B) The nonlinear region. Still isothermal in nature but now the viscous 
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element responds nonlinearly. At low De this portion of the traction curve 
can be described by a suitable nonlinear viscous function alone, while at high 
De a linear elastic element interacts with the nonlinear viscous element. 

(C) At yet higher values of slip the traction decreases with increasing 
slip and it is no longer possible to ignore the dissipative shearing and the 
heat that it generates in the film. Johnson and Cameron [1967] showed that 
the shear plane hypothesis advanced by Smith [1965] does account for most of 
their experimental observations in this region. More recently Conry et al 
[1979] have shown that a nonlinear viscous element together with a simple 
thermal correction can also describe this region. 


3-1 Isothermal Traction Analysis 


The rheological model that describes the traction under simple slip in 
all three regions of operation fairly well is the J & T traction model as 
presented by Johnson and Tevaarwerk [1977] : 

(3-1) | % * Hr) = 5 , 

where, 

r = shear stress, 

G = shear modulus, 

£ = time, and 
7= shear strain rate. 


[s] 

[1/s] 


The dissipative function F(t) is open to the choice of the researcher to fit 
the observed traction, but Johnson and Tevaarwerk [1977] found that the 
hyperbolic sine: 


(3-2) 



l 


where, 

r 8 = non-linear shear stress parameter. [Pa], 
tj - local viscosity [Pa.s], 


described all of their experimental results in regions (A) and (B) very well. 
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At higher pressures and for fluids with high traction coefficients this 
dissipative function may be replaced by the purely plastic behavior of the 
material : 

(3-3) F(r) =0 for t < t c ,* F(t) = 7 for r = t c , 


where, 

r c = limiting shear strength of the fluid [Pa]. 

Whether the perfectly plastic behavior of the material is intrinsic is not 
clear. Work by Johnson and Greenwood [1980] suggests that it is possibly the 
result of thermal behavior of the sinh model. For many applications the 
elastic/plastic traction model is adequate. It was used by Tevaarwerk and 
Johnson [1979c] and Tevaarwerk [1979b] to predict traction under various 
conditions of slip and spin. The analysis is completely isothermal in nature 
and for simple slip the traction is given by: 


(3-4) 



tan‘*S + — - — 

(1 + S 2 ) J 


I 


where, 


S = 


2 J1 

7 Tic 


J1 


{ ab 

“FT 


AU 

T 


and J4 = 


& 


The shear strain rate in the fluid was taken to be the same everywhere in 
the contact and assumed to be constant throughout the film thickness. Its 
magnitude was taken to be: 


7 = 


AU 

IT 


Equation (3-4) results from the integration of stresses caused by the 
shearing of an elastic element of pressure-independent average-shear-modulus G 
and the plastic stresses proportional to the local Hertzian pressure. The 
predicted traction from an elastic/plastic model compares very well with the 
experimentally observed values for combinations of slip and spin, provided 
that the spin or slip are not too large. Large slip or spin results in almost 
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purely dissipative stresses over the contact area and hence non-isothermal 
behavior. Traction prediction under these conditions is still possible but 
the thermal effects need to be brought into the picture. Tevaarwerk [1981a], 
[1979c] presents two techniques for calculating such spin traction curves. 

The former technique requires the shape of the traction curve in the large 
slip regime to provide a simple correction to the isothermal ly predicted spin 
traction. 

3-2 Thermal Traction Analysis 

The ability to separate the elastic stresses from the plastic ones can be 
used to perform a thermal traction calculation. The analysis presented here 
follows the technique outlined by Tevaarwerk [1980] and [1983]. 

Equation (3-4) resulted from the integration of isothermal elastic and 
plastic stresses over the contact area of an ellipse. For the region of 
contact under elastic stress, the shear energy is conserved and therefore does 
not give rise to temperature increases. The plastically deforming region, 
however, is non-conservative and a temperature rise in the fluid is expected. 
This will lead to a reduction in the local strength of the fluid. Equation 
(3-4) may therefore be better written in its elastic and plastic parts: 


(3-5) 04 = J4 p + J4 # 


(3-6) 

(3-7) 


J4. 


4 S 

*(1+S 2 ) 2 


, and 



[ tan^S 


, S(S 2 -1 )1 

(1+s 2 ) 2 J 


Equation (3-5) may now be modified by r e /r 0 where r e is the average 
stress under thermal conditions and r 0 is the average stress under isothermal 
conditions. We are dealing therefore with averaged stresses in the 
plastic region of the contact even though the isothermal stress distribution 
is according to the Hertzian pressure. This seemingly contradictory 
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assumption is supported by theoretical evidence by Tevaarwerk [1979c]. The 
modified equation would therefore be: 


(3-8) J4 t = J4 e + ^ j 4p 


The modification term r c /r 0 can be found from a thermal balance over 
the contact region under plastic stress. This region can be thought of as a 
thermal source whose heat is conducted/convected away. The length of the 
source is a function of the location of the onset of plastic deformation after 
the initial elastic region, however in this simple model we will take the 
source length to be "a" where this is the semi contact length in the running 
direction. As a simple thermal balance we will use the expression reported by 
Johnson and Cameron [1967] for the shear plane temperature: 


(3-9) 


where , 


*c - *0 - q £ 


F, + 


.5 ' 
JP e J 


t 


(3-10) 



[-] 


and, 


q = average thermal strength of the source, 

8 C = shearplane temperature in the contact, 

6 0 = inlet temperature of the fluid, 

k f = thermal conductivity of fluid in the contact, 

h = central film thickness in the contact, 

k s = thermal conductivity of the roller material, 

C 8 = specific heat of the roller material, 

p s = density of the roller material, and 

F s = thermal resistance of the film. 


'W/m 2 ] 
°C] 
r°C] 


’C] 


[W/m a 

■] 

W/m°( 
iJ/kg^C] 
: kc|/m 3 ] 


>C] 


This expression is valid only when the heat is conducted through the film and 
convected away by the discs, a condition that is true for most concentrated 


contacts. The strength q of the source is given by: 
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(3-11) q = rAU . 

The thermal resistance of the film F, may be found by calculating the average 
film temperatures due to distributed sources and is normally taken to be: 


(3-12) F, = 0.1. 


In order to proceed any further we need a relationship between temperature 
and the shear strength of the fluid. A typical relationship that has its 
roots in the Eyring theory of fluid transport is given by: 


(3-13) 


m = 


1 

r 


A + BP + (0+D) Ln 



I 


where, 


A = viscosity temperature constant, 

B = pressure viscosity constant, 

C = non-linear shear stress constant, 

D = fluid solidification temperature, 

E = fluid viscosity, and 
P = pressure of the fluid in the contact. 


°C] 
°C/Pa 
°C/Pa 
„°C] 
[Pa.s] 
"Pa] 


At first sight it seems that this equation has five disposable constants in 
it. However, two of these constants (A and D) may be obtained from the 
atmospheric viscosity temperature relationship, and one more (B) can be 
obtained from the Barus viscosity pressure relationship. The constant D is 
known as the solidification temperature, the temperature to which the fluid 
should be cooled to become solid like under atmospheric pressure. Only the 
constants C and E need to be determined experimentally from the traction 
results and this will be done in the next chapter. It should be noted here 
that the ultimate aim is to derive the fluid traction parameters such that 
they apply for all the experimental conditions reported here. 

By using equations (3-9), (3-10), (3-11), and (3-13) the average thermal 
shear stress can be obtained for a given set of conditions. In equation (3-8) 
we need the ratio of the average contact shear stress under thermal conditions 
to that under isothermal conditions. This is really the ratio of the shear 
stress given by equation (3-13) evaluated at the shearplane temperature 6 C 
(from equation 3-9) and the shear stress as evaluated at the inlet temperature 


conditions 6 0 . The other contact conditions remain the same for this ratio 
calculation. 

The method that is outlined above was used for the traction data analysis 
as reported by Tevaarwerk [1985]. For the analysis of the traction data as 
reported here, several improvements were made in the model and these changes 
are outlined in the next chapter. 
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4-0 EXTRACTION OF THE TRACTION PARAMETERS 

In order to use the traction model to calculate the level of traction in 
a given contact condition, the relevant traction parameters need to be 
extracted. In equation (3-1) we saw that both the elastic effects and the 
nonlinear viscous effect account for the behavior of the fluid. Elastic 
effects are typically felt under conditions of high contact pressure, low 
temperature, and high speed. 

4-1 Extraction of The Shear Modulus Parameters 

As discussed in chapter 3 the initial linear slope on a traction curve 
can be the result of either a viscous or an elastic response of the material 
in the contact to the strain implied. The parameter that determines the 
actual response is the dimensionless grouping of the relaxation time of the 
material in the contact and the transit time of this material through the 
contact. For a simple Maxwell type material this number is known as the 
Deborah number and is given by: 


(4-1) De - $ 


where, 

ri = viscosity of the fluid in the contact, 

U = transit velocity of fluid in the contact, 
a = semi contact length in rolling direction, and 
G = shear modulus of the fluid in the contact. 


Pa.s] 

m/s] 

m] 

Pa] 


Under the assumptions of a simple pressure distribution according to Hertz, a 
constant film thickness in the contact, and a constant shear modulus over the 
contact area, the initial small strain behavior for the material will be 
elastic if the Deborah number is larger than 10 and will be viscous in 
response if the Deborah number is less than 0.1. In between these two values 
the response is due to a mixed viscoelastic behavior. 

Many assumptions have been made before arriving at this point, however 
even with more complicated analyses where the pressure is allowed to influence 


21 


the viscous and elastic properties it is found that the transition points 
occur at about the same Deborah numbers. Also for the range of parameters 
normally encountered in traction contacts the Deborah number is such that the 
initial slope of the traction curve is almost always governed by the elastic 
properties of the contacting material. In the analysis of the data as 
performed here this assumption is implicit. 

4-1-1 Shear Modulus For Constant Properties 

When the initial linear response is completely elastic it is quite easy 
to calculate the value for the actual modulus that caused this slope. Under 
the assumptions of constant properties throughout the contact this modulus can 
be extracted from the initial slope using the following equation: 


(4-2) G - m [Pa] . 

This equation is applicable regardless of whether the slope is measured under 
longitudinal slip or under side slip. Using the known values of m, P 0 , and a, 
and a calculated value of h for film thickness, it is quite simple to 
calculate the fluid shear modulus 6. 

4-1-2 Shear Modulus With a Simple Compliance Correction 

One of the criticisms that is often raised at the above analysis is that 
it neglects the influence of the disk compliance on the measured slope. Disk 
compliance is the result of the elastic creepage of the disk material due to 
the tractive stresses on the surface. The traction response in the initial 
linear range is affected by this disk creepage in that it makes the slopes 
lower than if the discs were infinitely stiff. An exact correction of the 
modulus for the disc compliance is not possible at the moment. The analysis 
that is presented here is that due to Johnson and Roberts [1974]. 

If we let m' be the slope of the traction curve for dry rolling bodies, 
then from the addition of the compliances of the discs and the film a simple 
corrective term for the shear modulus may be derived as shown in 
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Equation (4-3). From this expression it is obvious that when the measured 
slope approaches the dry slope the corrected value for the shear modulus tends 
to infinity. 


<«-3> G <‘ G (n PTft) ' [Pa] 

where G c = simple compliance corrected modulus. 

The dry slope m' can be calculated from the expression given by Kalker [1967] 
as: 


(4-4) m' = jjp- , [-] 

where, 

G 8 = shear modulus of the disc material, and [Pa] 

Q = Kalker coefficient. [-] 

The value of the Kalker coefficient depends on the aspect ratio of the 
contact and the direction of slip. Below is a table which gives these values 
for the tests as reported here. 


Table 4-1. Kalker Coefficients Q for Sideslip 

aspect ratio 

coefficient 

k 

Q 

1 

.56 

2 

.70 

5 

.81 
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4-2 Extraction of the Large Strain Parameters 

The large slip region , that is the region beyond the traction peak, is 
almost exclusively governed by the dissipative element in the rheological 
equation. With the elastic effects becoming insignificant it is now quite 
easy to extract the governing parameters for this region. In essence what is 
required is a reverse analysis of the traction calculation normally used for 
the calculation of the traction curves. When the elastic response of the 
fluid is no longer dominant, the shear strain rate in the film is taken by the 
viscous element and we may express equation (3-1) as: 

(4-5) F(t) = 7 = y sinh ( ^-) 


When the argument r/r, is larger than 1 this may be further simplified to 
the following: 

(4-6) F (r) = 7 = ^ Exp ( ^ ) . 

This expression relates the stress to strain response for a small element 

O 

of fluid that is subjected to a shear strain rate, 7 , at a given pressure and 
a given temperature. Each of the parameters in this equation can be a 
function of temperature or pressure or both. In the analysis of traction it 
is normally sufficient to treat the fluid viscosity as a function of 
temperature and pressure and to make the non-linear stress parameter a 
function of temperature only. The following relationships will be used here: 

— for the viscosity we will combine the Vogels temperature viscosity and 
the Barus pressure viscosity relationships as follows: 

1 <*• P) ■ v. exp ( jfis; +^sr) . 


( 4 - 7 ) 


[Pa.s] 
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where, 

rj 0 = viscosity constant, 

A = temperature viscosity constant, 

D v = solidification constant for Vogels equation, 
B v = pressure viscosity constant, 

D = solidification constant for Barus equation, 

0 P = local temperature, and 

p = local pressure. 


Pa.s] 



Pa] 


— for the temperature dependence of r s on the local temperature we will 
use the following: 


C = non linear shear stress constant [°C/Pa] , and 
D, = solidification constant for Eyring equation [°C]. 

These parameters (except for the temperature and the pressure) are all 
constants for a given fluid. 

In order to proceed any further we have to make some assumptions about 
the conditions within the contact. First, we will assume that average 
conditions exist for the contact. Hence temperature and pressure are only 
mean temperatures and pressures. Second, we will assume that the important 
temperature for the analysis is the shear plane temperature as given by 
Equation (3-9). Third, the shear strain rate 7 will be assumed to be constant 
throughout the film and, in the case of the traction experiments, is constant 
over the contact area. Its magnitude is given by: 


(4-8) 

where , 


(4-9) 7 s Av/h, and [s ' 1 ] 

h = central filmthickness [m] 

The central film thickness is calculated from the expressions by Hamrock 
and Dowson [1977]. Further modifications to this film thickness to allow for 
inlet shear heating are made by using the Murch and Wilson [1975] approach. 
Hence in total we will make the following simplifications: 
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(4-10) 


local average equation # 


temperature 

e 


(3-9) 

shear strain rate 

0 

7 

Av/h 

(4-9) 

pressure 

P 

2P 0 /3 

(2-1) 

shear stress 

T 

F y /*rab 



By substituting the expressions from equation (4-7) to (4-9) into Equation 
(4-2) we obtain: 


(4-11) 


AV 0 C + D S r- ( T c C A BP ^ 

IT = — 2CV7 Exp V"57D # " *^+D v " "37R) p ') 


At a first glance there appear to be 7 parameters that can be used to fit the 
data. This is not the case however because several of these are determined by 
the results from other experiments. For example, the correlation of the 
viscosity data for the lubricants to equation (4-3) solves for the constants 
A, B, D v , and D p . It remains therefore to find the constants D s and C. These 
can be derived from the thermal region of the traction curve itself by curve 
fitting Equation (4-11) to it. For this purpose it is better to write this 
equation in a slightly different form: 


(4-12) 




BP 


+ Ln 


( 


2CV„ Av 

M.) h 



From the above equation it is apparent that if this relationship holds then 
the results from the traction measurements should form a straight line when 
plotted in the above fashion. The slope of this line gives the value for C 
while the value of D # can be calculated by adjusting it until the best fit is 
obtained. 
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4-3 Reduced Pressure Effects 

When a fluid film is present in the contact zone the pressure 
distribution is not strictly Hertzian but modified by the fluid. Due to the 
hydrodynamic action the pressure distribution is more peaky, and is spread 
over a broader area. This will lead to a reduction in the mean pressure in the 
contact. The reduced pressure may be calculated from the pressure ratio as 
given below: 

(4-13) P r = P [1 - 4 G e - 0 25 YI' 3 e<' 2 3 / k) ] , 

where, 

G, = Johnson's elasticity constant, 

P = mean Hertz contact pressure, 

YI = Inlet shear heating factor, 
k = aspect ratio (b/a), 

P r = reduced contact pressure, 

4-4 Transit Time Dependence 

The fluid in the contact is subjected to very high pressure gradients fo 
a very short time duration. It is therefore not likely that the viscosity in 
the contact has reached equilibrium conditions. How near equilibrium that it 
is will depend on several factors such as the time it is subjected to the 
pressure and the ease with which the molecules can rearrange themselves. 

A suitable dimensionless grouping that takes these various factors into 
account is given as: 


E-f 

Pa] 


(4-14) * = . 


where, 


$ ■ dimensionless time delay parameter, 

U * mean roller speed, 

Tj(6 0 ) = viscosity at inlet temperature, 

a = semi Hertz contact length, and 
P 0 = Hertz contact pressure. 


] 

m/s] 

Pa.s] 

M 

'Pa] 
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This parameter is used in conjunction with the viscosity constant in the 
following way: 


(4-15) V 0 = E^, [Pa.s] 


where, 

E 0 = equilibrium viscosity, and [Pa.s] 

E 2 = time delay exponent. [-] 


The single constant V 0 has now been replaced by this two constants equation 
so we have increased our degree of freedom by one for the system. 

4-5 Thermal Conductivity Temperature Dependence 

It may be expected that the thermal conductivity of the fluid is a function of 
the temperature. This has a direct influence on the thermal resistance of the 
fluid in the contact. In order to take this effect into account the thermal 
conductivity is fitted to the following equation: 

(4-16) k f = 3 — j^p- . [W/m °C] 

Since the shear plane temperature occurs in equation (4-12) and the thermal 
conductivity occurs in equation (3-9), the solution to both the thermal 
conductivity and the shear plane temperature must be done on an iterative 
basis. 

4-6 Determination of the Experimental Parameters 

We are now in a position to use the data obtained from the traction 
experiments and fit the rheological model to them. As stated earlier, the 
viscosity temperature-pressure data for the traction data reduction may be 
obtained from isobaric and isothermal viscosity measurements. Similarly the 
thermal conductivity data can be obtained from published information on RP1 in 
the Liquid Propellant Manual [1982]. The following sections will deal with 
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the determination of the actual numerical values of the various constants 
needed for the equations of state and the constants needed for the fluid 
rheology. 


4-6-1 Determination of Viscosity Temperature Pressure and 
Temperature Conductivity Parameters 

Not all of the constants in Equation (4-12) need to be determined from 
the experimental traction data. The temperature viscosity characterization 
may be used to solve for A and D v while the pressure viscosity effects can be 
used to solve for B and D p . 

The viscosity temperature data may be obtained from the Liquid Propellant 
Manual [1982]. Pressure viscosity data for RP1 was taken from the data 
generated by Bridgeman [1970] for kerosene. Using these two sources of 
information allows the determination of the constants in Equation (4-7). 


(4-7) 


where, 


V (0, P) = V o ex P 




A = 542.4, 

B = 544, 

D v = 143, 

D d * 25.6, and 
rj 0 = 6.33 x 10‘ 6 


°C] 

GPa' 1 ] 

C°c] 

: c] 

[Pa.s] 


[Pa.s] 


Thermal conductivity data for RP1 was taken from the Vargaftik [1983]. 
Fitting this data to Equation (4-16), 

(4-16) k f = j- ^' k ,i . [W/m °C] 

results in the following values for k 1 and k" 


k' = 61.6, and 
k" = 535.5 


W/m °C] 

;°c] 
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4-6-2 Experimental Shear Modulus from Traction Data 

Similar to the non-linear stress parameter r c , it is found that the fluid 
shear modulus as extracted from traction tests is influenced by the inlet 
temperature, by the contact pressure, and by the rolling speed. As a fairly 
simple relationship the following expression is used: 


(4-17) 


„ G, + G, P + G 3 U| 



where, 


® 1,3 

* experimental constants, 


P 

= contact pressure, 


U 

= rolling velocity, and 


D v 

= Vogels fluid solidification temperature. 


[Gpa] 


GPa 

m/s 

’•C] 


Using the initial slopes from the traction experiments and some regression 
analysis the following values were found for the constants G 13 . 


G, 1.0 

G 2 = 0.25 

G 3 = -0.005 


"GPa °C] 

°C] 

‘GPa °C s/m] 


To obtain this data, the high pressure, low temperature curves were mostly 
used in the calculations. 


4-6-3 Experimental Traction Constants 

1 

, By combining the Equations (4-8), (4-9) and (4-10) the large slip traction 

i 

data may be fitted to the following equation: 


(4-12) 


iSfc 


1 

Z 




BP 


+ Ln 


( 2CV. av ) ' 
V (fl’+D.) h ) J 


I 


where the shear plane temperature is calculated from equation (3-9) as: 
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(3-9) 




P is given by Equation (3-10) as: 


(3-10) 



and F s is: 


(3-12) F, « 0.1 


The constants that remain to be solved for at this point are: 

C, D # , E 0 , and E 2 . 

These can now be solved for by using a successive regression analysis on the 
large slip traction data. This was performed using the experimental traction 
data as reported in Chapter 2. A summary of the test data used and the 
calculated values for some of the intermediate parameters are shown in 
Table 4-2. The thermal properties for the test disc material was taken to be: 

k s = 15, [W/m °C] 

(4-18) C s = 500. [J/kg °C] 


The degree of fit obtained from the data on Equation (4-12) is shown in Figure 
4-1. From this regression analysis the following constants resulted: 


C = 2.75 x 10‘ 5 ' 

D s - 143, 

E. = 6.96 x 10" 5 , and 
E 2 = -0.057. 


■°C/Pa] 

;°c] 

Pa.s] 

[-] 


(4-19) 
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TABLE 4-2. TRACTION CURVE DATA ANALYSIS FOR RP1 


The following oil properties were used in the analysis: 
Viscosity temperature 
Viscosity Pressure 
Thermal Conductivity 
Eyring Ogress- 25. E 


A= 542.4 ’C D-143.0 ’C V0- .0000533 Pas 
B= 544.0 ’ C/GPa Dp= 25.5 ’C 
K'« 535.5 W/m’C K’ 51.7 ’C 
Peak Factor- .98 


Disc Thermal Conductivity 


Ks 1 


15 W/m’C Specific Heat-500 J/kg’C 




TEST INCLUDED IN THE MULTIPLE TRACTION DATA ANALYSIS 


Test t 

Po 

Uo 

To 

Fz 

kk 

a 

Pes 

Tshear 

H0 

Yi 

Plot 

( > 

( GPa ) 

<m/ 

5 > 

( ’C> 

(N) 

(.-) 

( mm ) 

(-> 

( ’C ) 

mic 

— 

Mark 

RP196M1TR 

1 

.01 

10. 

3 

34 

200 

1.7 

.235 

2 

47 

.097 

.98 

A 

RP196M2TR 

1 

.01 

30. 

4 

35 

200 

1.7 

.235 

26 

61 

.193 

.92 

B 

RP195M3TR 

1 

.01 

49. 

6 

39 

200 

1.7 

.235 

65 

71 

.237 

.87 

C 

RP19EM4TR 

1 

.27 

11. 

5 

37 

400 

1.7 

.297 

2 

61 

.093 

.98 

D 

RP196M5TR 

1 

.27 

30. 

8 

38 

400 

1.7 

.297 

18 

76 

.175 

.92 

E 

RP196M6TR 

1 

.27 

50. 

0 

41 

400 

1.7 

.297 

44 

92 

.211 

.87 

F 

RP196M7TR 

1 

. B0 

11. 

9 

39 

800 

1.7 

.374 

1 

75 

.086 

.98 

G 

RP19BM8TR 

1 

.60 

30. 

5 

40 

800 

1.7 

.374 

11 

99 

.150 

.93 

H 

RP195M9TR 

1 

.60 

50. 

5 

44 

800 

1.7 

.374 

31 

117 

.188 

.88 

I 

RP19BMATR 

1 

.92 

11. 

8 

40 

1380 

1 .7 

.448 

1 

91 

.077 

.98 

J 

RP196MBTR 

1 

.92 

31. 

4 

44 

1380 

1.7 

.448 

7 

119 

.128 

.94 

K 

RP19BMCTR 

1 

.92 

50. 

0 

49 

1380 

1.7 

.448 

19 

141 

.157 

.89 

L 

RP196N1TR 

1 

.01 

49. 

5 

58 

200 

1.7 

.235 

37 

84 

.175 

.90 

M 

RP196N2TR 

i 

.01 

30. 

7 

57 

200 

1.7 

.235 

13 

74 

.135 

.94 

N 

RP196N5TR 

l 

.27 

30. 

9 

60 

400 

1.7 

.297 

9 

87 

.123 

.94 

0 

RP196N6TR 

l 

.27 

50. 

0 

62 

400 

1.7 

.297 

26 

97 

.158 

.90 

P 

RP196N7TR 

l 

. B0 

50. 

5 

54 

800 

1.7 

.374 

19 

122 

.146 

.90 

Q 

RP19BN8TR 

l 

.60 

31. 

9 

64 

800 

1.7 

.374 

6 

108 

.110 

.94 

R 

RP196N9TR 

l 

.60 

11. 

3 

64 

800 

1.7 

.374 

0 

91 

.054 

.99 

S 

RP196NATR 

i 

.92 

11. 

2 

64 

1380 

1.7 

.448 

1 

107 

.057 

.98 

T 

RP196NBTR 

l 

.92 

31. 

1 

56 

1380 

1.7 

.448 

4 

127 

.097 

.95 

U 

RP196NCTR 

l 

.92 

49. 

6 

69 

1380 

1.7 

.448 

12 

142 

.127 

.91 

V 

RP196NETR 

i 

.01 

30. 

5 

64 

200 

1.7 

.235 

11 

79 

.122 

.94 

w 

RP196NGTR 

i 

.27 

51. 

0 

63 

400 

1.7 

.297 

26 

99 

.159 

.90 

X 

RP196NHTR 

l 

.27 

30. 

.6 

51 

400 

1.7 

.297 

9 

89 

.121 

.94 

Y 




FIGURE 4-1. DEGREE OF FIT OF RP1 TRACTION DATA TO EQUATION (4-12) 
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5-0 TRACTION PREDICTION 

The ultimate aim in obtaining the fluid rheological properties is to 
solve for the traction response of a concentrated contact under general 
conditions of slip, spin, and geometry. In this report we will restrict 
ourselves to the conditions of sideslip or longitudal slip only. Calculation 
methods for other contact kinematic conditions are outlined in Tevaarwerk 
[1979a] and Tevaarwerk [1979b] . 

As another way to observe the faithfulness of the traction modeling it 
may be instructive to predict the measured traction curves on the reduced data 
as extracted. 

5-1 Comparison Between Measured and Predicted Traction 

The process essentially involves the reverse operation to that used for 
the extraction of the traction parameters. Starting with Equation (3-8) 


(3-8) J4 t = J4 e + i J4 p 


where, J4 e = 


4 S 

*(1+S 2 ) 2 ' 


J4 


p 



tan _1 S + 


S(S 2 -1) I 

(1+S 2 ) 2 J ' 


and 


$ = 2 J1 = 2 G a AU 

~ 3 Ik ~ 3 t 0 h U 

As an average shear modulus G we have to use the combination of fluid 
shear modulus and the elastic compliance of the discs. This apparent modulus 
can be calculated from the compliance corrected modulus as given in 
Equation (4-17) together with the relationships (4-2), (4-3), and (4-4). This 
results in: 
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(5-1) 


G = 


1 + 


G c 

8Qa G c 

~imr s 


[Pa] 


where. 


G c 

a 


G 

Q 


s 


h 


compliance corrected modulus as given by Equation (4-17) 
semi contact size in rolling direction, 
shear modulus of disc material, 

Kalker coefficient, and 
central contact film thickness. 



The average fluid shear strength r 0 based on the inlet temperature can be 
calculated directly from Equation (4-12) by using 0 O as the shear plane 
temperature, 


(5-2) 




BP 


+ Ln 


( 


2CV„ Av 
(0 O +D S ) h 


) 


and 


using the fluid constants as given in Equation (4-19). 

To solve for the ratio of fluid shear strength at the shear plane temperature 
0 c and the inlet temperature 0 O requires a simple iteration since 0 C is a 
function of t c and vice-versa as indicated by Equations (3-9), (3-10), (3-11) 
and (4-12). Additional to the temperature dependence of r c through 0 C there 
is also the influence of thermal conductivity. In most cases only two or 
three iterations are needed for convergence. By using this technique, and the 
rheological data as extracted, the sideslip traction can be predicted for each 
of the measured traction curves. The comparison between measured and 
predicted traction may be observed from the graphs included in Appendix II. 
Overall the degree of predictions correlates very well with the experimentally 
observed results indicating a high degree of confidence in the traction 
model i ng . 
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5-2 Traction Predictions for Line Contact 

Traction prediction for line contact is essentially the same as that 
outlined in Section 5-1. The fundamental fluid parameters are independent of 
the exact contact geometry so there would be no change there. Changes do need 
to be made in the following: 

1) Film thickness calculations, 

2) Kalker correction factor, 

3) Mean contact pressure, and 

4) Shear plane temperature calculation. 

The film thickness calculations should reflect the increased film that is 
found when rollers form a line contact. In the line contact calculations 
performed here, the central film thickness calculations according to Dowson- 
Higginson [1977] are used. The effect of the increased film thickness is a 
slight reduction in the traction in the thermal region because of increased 
thermal resistance of the thicker film as compared with elliptical contacts. 

Kalker' s correction factor for the compliance of the discs should now 
reflect the direction of the slip and the contact geometry. From Kalker 
[1967] we find that a value of Q for line contact is: 

(5-3) Q = 0.6. [-] 

When the traction data was fitted to the fluid rheological model, averaged 
quantities were used. The average pressure in a line contact is somewhat 
higher for a given peak contact pressure than in an elliptical contact. From 
Hertz stress considerations we can calculate a correction factor for line 
contact. Hence for the same peak Hertz pressure we should use a mean pressure 
that is about 18 percent higher in the traction calculations. 

The shear plane calculations are based on Equation (3-9). This equation 
consists of two thermal resistance terms, 1) that due to the film, and 2) that 
due to the contact. The change in film resistance is taken care of by using 
the film thickness calculations for line contact. The change in contact 
resistance can be corrected for by using the ratio of contact resistance for 
circular contact and rectangular contact. Jaeger [1942] gives this ratio as: 
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-£li« = 1.05. 

''circular 

This is based on average contact temperature. Using the corrected terms in 
relevant equations, calculations identical to these for elliptical contact 
under side slip can be carried out. 

5-2-1 Line Contact Traction Predictions 

One of the aims of this investigation was to provide traction data on RP1 that 
can be used in the bearing analysis code SHABERTH. The traction data required 
by SHABERTH is in the form of a number of selected points from line contact 
traction test. The above outlined technique was used to modify the data 
obtained using elliptical contact traction tests into line contact data. 

Using this method the following line contact traction test curves were 
predicted. 


Contact Pressure 

P 0 : 

1.6, 2.3, 3.0, 

[GPa] 

Rolling Speed 

U : 

30, 60, 100, 

m/s] 

Temperature 


40, 75, 100, and 

;°c] 

Slide/Roll Ratio 

au. 

TP 

0.0 to 0.40, 1 

[-] 


The predicted line contact traction curves are included in Appendix III. 
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ORIGINAL MEASURED TRACTION DATA 
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APPENDIX II 


COMPARISON BETWEEN MEASURED AND 
PREDICTED ELLIPTICAL CONTACT DATA 




Po- 1.01 GPa s Tin- 35 'C : U- 10 m/s 







Traction Coefficient Fy/Fz l-l Traction Coefficient Fy/Fz 



Sidesl ip DV/U C-D 



Sideslip DV/U C-3 




Pred: # f Po a 1,27 GPa 











efficient Fy/Fz 













Traction Coefficient Fy/Fz [-] Traction Coefficient Fy/Fz 




Sides 1 ip DV/U [-] 






APPENDIX III 


PREDICTED LINE CONTACT DATA 




Traction Coefficient Fx/Fz C-] Traction Coefficient Fx/Fz 



Slide/Roll Ratio DU/U Z-l 



Slide/Roll Ratio DU/U C-D 



Traction Coefficient Fx/Fz C-3 Traction Coefficient Fx/Fz 




Traction Coe-f f 1 c 1 ent Fx/Fz [-] Traction Coefficient Fx/Fz 
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Pred : 
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Pred: 
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,Po- 3.00 
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U- 

30 
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Pred : 

+ 

# Po“ 3.00 
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U- 

60 

m/s 

Pred : 

X 

f Po“ 3.00 
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U- 100 

m/s 
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